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Abstract. - 

We consider a gas of neutral fermionic atoms at ultra-low temperatures, with the attractive 
interaction tuned to Feshbach resonance. We calculate, the variation of the chemical potential 
and the energy per particle as a function of temperature by assuming the system to be an ideal 
gas obeying the Haldane- Wu fractional exclusion statistics. Our results for the untrapped gas 
compare favourably with the recently published Monte Carlo calculations of two groups. For 
a harmonically trapped gas, the results agree with experiment, and also with other published 
work. 



Consider a dilute gas of neutral fermionic atoms at a low temperature. In general, the 
low-energy properties of the gas are determined by the scattering length a, the number density 
n, and the temperature T of the gas ( the effective range rg is small, so that ro/|a| ^ as a 
becomes large ) . When the attractive interaction between the atoms is increased continuously 
by magnetic tuning from weak to strong, the scattering length a goes from a small negative 
to a small positive value. In between, there is a zero-energy two-body bound state, and \a\ 
is infinite. The gas is said to be at unitarity in this situation, and the length scale a drops 
out. The behaviour of the gas is expected to be universal at unitarity [1]. Experimentally, if 
the temperature is small enough, a BCS superfluid is observed at the weak end, and a EEC 
condensate of dimers at the strong end [2] . This was predicted long back by Leggett [3] , who 
extended the BCS formalism in a novel fashion to analyse the physical situation. The BCS to 
BEC transition is found to be smooth, with no discontinuity in properties across the unitary 
point. There has been much interest amongst theorists to calculate the properties of the gas 
in the unitary regime (fcy|a| >> 1), where kf — (Stt^ti)^/'^ is the Fermi wave number of the 
noninteracting gas. This is a challenging task, since there is no small expansion parameter, 
and a perturbative calculation cannot be done. In particular, at T = 0, the energy per particle 
of the gas is calculated to be 

^ (1) 
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where ^ ~ 0.44 [4]. The experimental value is about 0.5, but with large error bars [5]. Recently, 
there have been two Monte Carlo (MC) finite temperature calculations [6,7] of an untrapped 
gas at unitarity, where various thermodynamic properties as a function of temperature have 
been computed. For a harmonically trapped gas, there are experimental results [8], as well as 
theoretical calculations [9]. 

In the unitary regime, the thermodynamic properties have both bosonic and fermionic 
features [6], and it is natural to ask if in this regime the quasi-particles obey a statistics 
which is intermediate between the two. In this paper, we suggest on general grounds that 
at unitarity, so far as average properties of the system are concerned, it should behave like 
an ideal gas obeying the generalised exclusion statistics of Haldane [10]. The definition of 
the statistical parameter, denoted by g{> 0) in the present paper, is based on the rate at 
which the number of available states in a system of fixed size deceases as more and more 
particles are added to it. The statistical parameter g assumes values and 1 for bosons and 
fermions respectively, because the addition of one particle reduces the number of available 
states by g. We first deduce the value of g for the unitary gas from theory using Eq.(|T]l. 
fitting ^ = 0.44. The value of g thus determined remains the same independent of the nature 
of confinement as it should since the microscopic origin of the value of g depends only on 
the interaction between fermions and not on how the system is prepared. The application of 
the finite temperature distribution function [11] (called Haldane- Wu statistics in this paper) 
then enables us to calculate the temperature dependence of the energy per particle, and the 
chemical potential of the unitary gas. Our results for both trapped and untrapped gases are 
in good agreement with experiment, and MC calculations. 

We shall now give the rationale for using Haldane- Wu statistics at unitarity: 

• A strong hint that this may be the case comes from the observation that the kinetic and 
potential energies scale the same way when there is no length scale left from the inter- 
action. As is well known, Haldane- Wu statistics is realised by the Calogero-Sutherland 
model in one dimension [13]. The potential and kinetic energy both scale the same way 
in this model, and both the energy densities scale as n^. Similarly, fermions in two di- 
mensions interacting with a zero-range potential have their kinetic and potential energy 
densities scale as n^, obeying Haldane- Wu statistics [14]. 

• In the present case, a compelling evidence comes from the fact that the second virial 
coefficient of the gas at unitarity is temperature independent [15]. In exclusion statistics, 
the scale-invariant interaction between atoms alters the ideal Fermi (Bose) values of the 
(exchange) second virial coefficient +{—)2~^/^ by adding an interacting part [16]. 

The above arguments are heuristic and indicative. A quantitative understanding can be 
obtained only when the effective interaction is known fully. In the absence of such a theory, in 
this paper we pursue a phenomenological approach where we assume the validity of exclusion 
statistics on the average for quasi-particles which are otherwise non-interacting. The effect 
of interaction is entirely subsumed in defining the statistics of the quasi-particles. We first 
estimate the value of the statistical parameter g from the following considerations: 

For Haldane-Wu statistics, the distribution function (or occupancy factor) in a single 
particle state with energy a is given by fi = {wi + g)^^ , where Wi obeys the relation 

w!il + w,y-3 ^exp[ie,~fi)p] , (2) 

where f3 = 1/T, T being the temperature in units of the Boltzmann constant. Note from the 
above that for g = and 1, the distribution function fi reduces to the familiar bosonic and 
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fermionic forms. It is also clear that for T = 0, the occupancy factor is 

/.(T = 0) = - , e, <tJL, 
9 

MT = 0) = , e^> II . 



(3) 



Now consider N fermionic atoms obeying this statistics at T = in a large volume V . The 
new Fermi momentum kf is determined from the relation 



9 (27r)3 



where we have included a spin degeneracy factor of 2. The modified Fermi momentum kf, 
from above, is fc/ = g^^^kf, where kj is the fermi momentum of the noninteracting Fermi gas. 
It also follows that the energy per particle of the unitary gas is given by 

E ^ /3 3^ 
iV 5 2M ■ 

Comparing with Eq.l^, we see that ^ = 5^^"^, and choosing g — 0.29 gives the generally 
accepted value of ^ = 0.44. This therefore fixes the only free parameter in the model, namely, 
g and it should be valid independent of temperature and the nature of confinement as it is 
the parameter which determines the statistics of quasi-particles. 

The main advantage of our model, however, is the calculation of the bulk properties of 
the gas as a function of the temperature, and this we proceed to do now. We follow the well 
known method (see for example the paper by Aoyama [17]) for this purpose. For a given 
density of single-particle states -D(e), we have 



N = 



D{e)de 

[w + gY 



E = 



eD{e)de 

{w + g)' 



For the 3— dimensional gas, D{e) — C-^/e, where the constant C = |A^e^^^^. 



(4) 

Furthermore, 



Cf = r^j^j is the Fermi energy of the noninteracting Fermi gas. Changing the variable from 
dt to dw, and using the relation involving lu's given above, one gets after some algebra 




(5) 



(6) 



In the above, wq is the value of ui at e = 0. For our choice oi g = 0.29, the Eq. ^ is solved 
at a given (T/ef) for wq numerically, and this wq is used in Eq.@ next to obtain [E/Ntf). 
From the definition of , it also follows that the chemical potential /i at temperature T obeys 
the relation 



— = [g In Wo + (1 - g) ln(l + wo) 



(7) 



Our results for the energy per particle and the chemical potential (in units of the noninteracting 
Fermi energy tf ) are plotted in Fig^and Fig|21 respectively. 
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Fig. 1 - Plot of the energy per particle as a function of temperature. Both the abscissa and ordinate 
are in units of the free Fermi energy. The solid line corresponds to our calculations with g — 0.29. 
The solid squares (green) with error bars are the MC calculations of [7], and the triangles (blue) are 
the MC calculations of [6]. 



Our results are not sensitive to the fine-tuning of the statistical parameter g. In Fig. ^ we 
also show, for comparison, the recent MC calculated points of Bulgac et al. [6] and Burovski 
et al. [7]. It will be seen that the agreement is very good, although the chemical potential n 
as calculated by us starts to differ from Burovski et al result for T/e/ > 0.8). 

The finite temperature results are easily generalised for fermions in harmonic trap. Con- 
sider the fermions at T = 0. The density of states D{e), including a spin degeneracy factor 
of 2, is /(TiLoY, where the oscillator parameter is defined as w = {ujxUiyUJzY^^ ■ It follows 
immediately that e/ = {"igNY/^fiuj^ and the energy E = g^/^{?>N)'^/^ / A Tito. These results are 
the same as the Thomas- Fermi density functional approach of Papenbrock [18]. We can easily 
extend these results to finite temperatures using this density of states in Eq.l^J. 





/•°° dw 


ln| 




J wo w{l+w) 







E 



dw 



w{l + w) 




(8) 



(9) 



The expression for ^ remains the same as Eq.(|7|l, although the numerical values of wq as a 
function of T are quite different from the unconfined gas. We present the results for average 
energy in the trap in Fig|31 using the same value of 17 = 0.29 since the statistical parameter 
depends only on the mutual interaction and not on the nature of confinement. It will be 
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Fig. 2 - Same as in Fig.Q for chemical potential plotted as a function of temperature. 

seen that the agreement with experimental data of Kinast et al [8] as well as the many body 
calculation of Hu et al [9] is very good. 

Other thermodynamic quantities could be readily calculated. However we note that the 
model cannot yield the two- or many-particle correlation functions. In this regard, the situa- 
tion is similar to the one-dimensional Calogero-Sutherland model [12], which can be mapped 
on to a system of quasi-particles which obey Haldane- Wu statistics [13]. But this does not 
help in obtaining the correlation functions, for which the full many-body calculation has to 
be done. Moreover, the ideal Haldane- Wu gas cannot describe super-fluidity. Therefore, the 
main usefulness of the present approach is its ability to calculate the temperature-dependence 
of various bulk properties of a unitary gas with just one free parameter, namely the statistical 
parameter g. 

* * * 

We are extremely grateful to Evgcni Burovski, Joaquin Drut and Hui Hu for sharing the 
results of their calculations with us. We would like to thank Akira Suzuki for useful discussions, 
and for invaluable help in preparing the manuscript. One of the authors (R.K.B.) benefited 
greatly from many discussions with Duncan O'Dell. Thanks are also due to Jules Carbotte. 
This research was supported by NSERC of Canada. 

REFERENCES 



[1] Baker G.A., Phys. Rev. C60, 054311 (1999); Heiselberg H., Phys. Rev. A63, 043606 (2001); Ho 
T.-L., Phys. Rev. Lett. 92, 090402 (2004). 



6 



EUROPHYSICS LETTERS 



2.5 




2 



1.5 



LU 

lD 



1 



0.5 - 



FESg=0.29 

Kinast etal ■ 
Huetal 











0.1 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 



T/Tp 



Fig. 3 - Plot of the energy per particle as a function of temperature for the confined gas. Both 
the abscissa and ordinate are in units of the free Fermi energy. The solid line corresponds to our 
calculations with g — 0.29 and the dashed lines corresponds to the calculations presented in Ref. [9]. 
The experimental data is taken from Kinast et al [8]. 
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